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ABSTRACT 

The millimeter- wave rotational transition line ( Jk = lo ~ Oo) of deuterated 
methane CH 3 D has tentatively been detected toward the low-mass Class pro- 
tostar IRAS 04368+2557 in L1527 with the Heinrich Hertz Submillimeter Tele- 
scope. This is the first detection of CH3D in interstellar clouds, if confirmed. 
The column density and fractional abundance of CH3D are determined to be 
(9.1 ±3.4) x 10 15 cm -2 and (3.0 ± 1.1) x 10~ 7 , respectively, where we assume the 
rotational temperature of 25 K. The column density and fractional abundance of 
the gaseous CH 4 are estimated to be (1.3-4.6) x 10 17 cm" 2 and (4.3-15.2) x 10~ 6 , 
respectively, by adopting the molecular D/H ratios of 2-7 % reported for various 
molecules in L1527. The fractional abundance of CH 4 is higher than or compa- 
rable to that found in high-mass star-forming cores by infrared observations. It 
is sufficiently high for triggering the efficient production of various carbon-chain 
molecules in a lukewarm region near the protostar, which supports the picture 
of the warm carbon-chain chemistry. 

Subject headings: ISM: individual (L1527) ISM: molecules 

1. Introduction 

Methane (CH 4 ) is the most fundamental hydrocarbon molecule in space, and is thought 
to play an important role in interstellar chemistry and planetary chemistry. The gas phase 
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CH4 has been identified toward bright infrared sources like high-mass young stellar objects 
by observing the infrared vibrati on-rotation spectra of the (C-H stretching: 3 lira) an d 
u 4 (CH 2 bending: 7 fan) bands flLacv et al. 199l[ iBoogert et al. 2004 knez et al. 2009h . 
Such observations are also carried out toward comets, planets, and satellites like Titan (e.g. 
Penteado et al. 20051 ). Furthermore, CH4 is known as one of the major constituents of icy 
grain mantles in interstellar clouds. The infrared absorption of the u% (3.3 /iin) and 
(7.7 urn) bands have been observed toward various sources including low- mass protostars 
(e.g. lLacy et al. 1991c lOberg et al. 20081 : lOberg et al. 20 111 ). In spite of this progress, our 
understanding of the distribution and the abundance of CH 4 in the gas phase is still limited 
particularly for cold interstellar medium including low-mass star-forming regions. This is 
because CH4 has no permanent dipole moment, and does not emit the pure rotational spec- 
tral lines which can readily be observed with radio telescopes. Furthermore, the low-mass 
protostars are not bright enough for high-resolution infrared absorption spectroscopy. 

One practical way to overcome this situation is to observe its mono-deuterated species, 
CH3D, which emits the pure rotationa l spectral lines due to its small electric dipole mo- 
ment (0.005657 D) JWofcv et al. 197p[ IWatson et al. 1979h . However, the past efforts to 
dete ct the Jk = In ~ On line (232.6 GHz) of CH3D were not successful even toward Orion- 



KL fjPickett et al. 1980c IWomack et al. 19961 ). Since active high-mass star-forming regions 
contain many kinds of saturated complex organic molecules, their contamination make it dif- 
ficult to detect the line. At that time, it was thought that detection of the CH 3 D line would 
be difficult in low-mass star-forming regions because of its relatively low column density in 
comparison with high-mass star-forming regions. 

Recently, it was recognized that CH4 plays a crucial role in chemical processes occur- 
ring in low-mass star forming regions. Sakai et al. (2008; 2009a) found low-mass star- 
forming regions which harbor extremely rich carbon-chain molecules. They are L1527 in 
Taurus and IRAS 15398-3359 in Lupus. According to the interferometric observation to- 
ward L1527, carbon-chain molecule s and their relat ed species show a steep increase in 
abundance inward of 500-1000 AU (jSakai et al. 20101), where the gas kinetic temperature 
is higher than 20 K (jShirley et al. 2002c Ijprgensen et al. 20021 ). Furthermore, they reside 
even in the gas infalling to the protostar. In order to explain the observational results, 
it is proposed that carbon-chain molecules are regenerated near the protostar triggered by 
evaporation of th e solid CH4 on dust grains (Warm Carbon Chain Chemistry: WCCC ) 
( jSakai et al. 20081 ). The desorption temperature of CH4 is about 25 K (jCollings et al. 2004 ). 
which is lower than that of H 2 (~ 100 K). Hence CH 4 can be evaporated in a luke- 
warm region (20 — 30 K) around a newly born protostar. The CH 4 evaporated from dust 
grains reacts with C + to form C2H^, which gives C2H2 and C 2 H by electron recombi- 
nation reactions. Further reactions of C2H2 and C2H with C + produce longer carbon- 
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chain molecules. A basi c part of this scheme is confirmed by chem ical model calcula- 



tions flAikawa et al. 20081 : lHassel et al. 20081 : lHarada and Herbst 20081 ). The WCCC pic- 



ture is consistent with the observed distribution of the carbon-chain molecules in L1527 
f lSakai et al. 201ph . 

Because of a high expected abundance of CH 4 , the WCCC sources would be a good tar- 
get to search for the CH 3 D line. The molecular D/H ratios are measured to be consistent for 
various molecules toward L1527 ( jSakai et al. 2009bl ). and hence, we can make a reasonable 
estimation for the abundance of gaseous CH 4 from the observation of CH 3 D. The result will 
verify the WCCC mechanism, and will lead us to a deeper understanding of a role of CH4 
in interstellar chemistry. With this motivation, we conducted a sensitive observation of the 
CH 3 D line toward L1527. 



2. Observations 

We first observed the J K = 1 - line (232.644301 GHz) of CH 3 D toward the low- 
mass protostar, IRAS 04368+2557 in L1527 with the IRAM 30 m telescope in April 2008. 
The observed position was (a 2 ooo, £2000) = (04 h 39 m 53 s .87, 26°03'09".7). We employed the 
A230/B230 receiver as a frontend. The beam size is 10. "6. In this observation, we found a 
hint of the CH 3 D line at the right Vlsr velocity with the intensity of 20 mK (Tmb), although 
the confidence level is only 2.3 a. 

Encouraged by this result, we then conducted a long integration observation of this line 
toward the same position with Heinrich Hertz Submillimeter Telescope (HHT) in January 
2011 and February 2012. The ALMA Band 6 (1.3 — 1.1) mm prototype receiver was used 
as a frontend in the 4 IF dual polarization mode, whose system noise temperatures ranged 
from 180 K to 250 K. In the observation in 2011, the CH 3 D J# = 1 — Oo line was placed in 
the lower sideband with the CS J = 5 — 4 line (244.9355565 GHz) in the upper sideband. 
The sideband rejection was about -13 dB. The beam size and the main beam efficiency of 
the telescope are listed in Table [TJ In the observation in 2012, the CH 3 D Jk = lo — Oo 
line was placed in the upper sideband with the C 18 J = 2 — 1 line (219.5603541 GHz) in 
the lower sideband. The telescope pointing was checked every 1.5 hours by observing the 
continuum emission of Jupiter, or the CS( J = 5 — 4) emission of S231 and IRC+10216, and 
was maintained to be better than 10". A small daily variation of the intensity was calibrated 



1 The HHT is operated by the Arizona Radio Observatory (ARO), Steward Observatory, University of 
Arizona. 
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by using the CS or C 18 line. The backend used was filterbanks whose bandwidth and res- 
olution are 64 MHz and 250 kHz, respectively. The frequency resolution corresponds to the 
velocity resolution of 0.32 km s -1 , which is comparable to or slightly smaller than the line 
width in this source (0.3-0.6 km s" 1 : Sakai et al. 2008; 2009). The observation was made in 
position-switching mode with the off-position of Aa = 20 arcminutes. The final spectrum 
was placed on the Tmb scale using observations of Jupiter. 



3. Results 

In the 2011 observation with HHT, we tentatively detected the Jk = lo ~~ Oo line of 
CH3D toward the low-mass protostar IRAS 04368+2557 in L1527, as shown in Figured! The 
line can be recognized in the both (H and V) polarization data, although it is marginally 
seen in the V polarization data due to a larger noise level. The line is not well resolved due 
to a limited frequency resolution of the backend employed. In preparation of the spectra, 
each 5 minute integration scan were carefully checked to eliminate the bad scan data with a 
heavy baseline distortion and/or an apparent strong noise spikes, as a usual data reduction 
process. Note that no apparent spikes appear at the line channel. For the total spectrum, 
the confidence level of the detection is 7.9 a in the integrated intensity, and 5.2 a in the 
peak intensity. This line shows the strongest intensity in the observed band (128 spectral 
channels), and the noise reveals a Gaussian distribution, as shown in Figure [2j Hence, the 
line is significantly detected in the total spectrum. The noise distributions of the H and 
V polarization spectra are also shown in Figures [2b and [2b, respectively. They show the 
Gaussian form, although the V polarization data have a slightly poor shape. 

To confirm the line, we also carried out the observation with the different frequency 
setting in 2012, where the CH3D line is set in the upper sideband, and marginally confirmed 
the line with the 3.6 a confidence level. However, we unfortunately suffered from serious 
spurious signals probably due to instability of the receiver and the backend in the 2012 
observation in contrast to the 2011 observation, and hence, we did not use the 2012 data for 
the final spectrum. 

The linewidth and the LSR velocity are consistent with those of other carbon-chain 
molecules observed in this source (e.g. Sakai et al. 2008; 2009b). We car efully checked the 



molecular line databases such as CDMS (jMuller et al. 20051 ). JPL Catalog (IPickett et al. 19981 ). 



and confirmed that no appropriate line exists at the frequency of ±0.5 km s -1 except for 
the CH3D line. We also confirmed no strong line in the image band. The lines of saturated 
complex organic molecules like HCOOCH 3 , C 2 H 5 CN, and (CH 3 ) 2 0, which give congested 
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spectra, are not detected in L1527 even with a very sensitive observation (ISakai et al. 20071 ). 
This is in contrast to the high mass star- forming regions like Orion KL, and also to the low- 
mass star-forming regions with the hot corino activity (e.g. IRAS16293— 2422; Cazaux et 
al. 2003). Furthermore, the linewidth in L1527 is narrower than those in Orion KL (~5 km 
s _1 ) and IRAS16293— 2422 (~2 km s _1 ). Although L1527 harbors long various carbon-chain 
molecules and their isomers, their rotational transitions in the 200 GHz region generally 
have very high upper state energies (>100 cm -1 , for carbon-chains with for or more carbon 
atoms) and are difficult to be excited in L1527. Therefore, an accidental matching of other 
lines would be unlikely, and our identification of this line to CH3D is the most reasonable. 

The column density of CH3D is estimated to be (9.1 ± 3.4) x 10 15 cm -2 from the 
integrated intensity by assuming the LTE (local thermodynamic equilibrium) condition with 
the excitation temperature of 25 K, which corresponds to the evaporation tem perature of 



CH4 from dust grains. The source size is assumed to be 20" (ISakai et al. 20101 ). Here the 
error is calculated from three times the standard deviation of the integrated intensities. The 
fractional abundance of CH3D is evaluat ed to be (3.0 ± 1.1) x 10~ 7 , where the H2 column 



density is assumed to be 3 x 10 22 cm 2 (IJorgensen et al. 20021 ). Temperature dependences 



of the column density and the fractional abundance are carefully inspected, as summarized 
in Table |2j The column density and the fractional abundance vary roughly linearly with the 
temperature over the range from 15 to 35 K. In this temperature range, the difference is 
within a factor of 2. 

varies within a factor of 2. The CH 3 D column density is much less than the upper limit 



reported for Orion-KL (N < 6 x 10 18 cm 2 ) by three orders of magnitude (IWomack et al. 19961 ) 



When we assume the deuterium fractionation ratio of 2-7 % by referring the ratios for var 



ious carbon-chain molecules observed toward L1527 (ISakai et al. 2009bl ). the CH 4 column 
density is evaluated to be (1.3-4.6) xlO cm" 2 by using the CH3D column density for the 
excitation temperature of 25 K. Hence, the fractional abundance of CH 4 is (4.3-15.2) xl0~ 6 . 
In this calculation, we do not explicitly consider the statistical factor due to four equivalent 
hydrogen atoms in CH4, because the deuterium fractionation ratios observed in L1527 are 
within the above range regardless of the number of equivalent hydrogen atoms (i.e. 0.05, 
0.03, 0.07, and 0.04 for 1-C 3 D/1-C 3 H, C 4 HD/C 4 H 2 , c-C 3 HD/c-C 3 H 2 , and NH 2 D/NH 3 , re- 
spectively). 
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4. Discussions 

In this study, we tentatively detected the rotational spectral line of CH3D in interstellar 
space for the first time. Although no cataloged lines do not match to the observed line, 
it might be a line of unknown species. To rule out this possibility completely, we need 
to observe the other rotational transitions. The intensities of the Jk — 2o — lo (465 GHz; 
E u = 33 K) and 3q — 2q (698 GHz; E u = 67 K) lines are predicted to be 21 mK and 12 mK by 
assuming the excitation temperature of 25 K and the linewidth of 0.5 km s -1 . However, their 
sensitive observations are difficult because of the relatively heavy atmospheric absorption at 
these frequencies. We have therefore decided to publish the present result as the tentative 
detection of CH 3 D, considering its importance in astrochemistry. 

The gas-phase CH4 abundance has been reported for a few high-mass star-forming 
regions by observing the infrared absorption lines. Lacy et al. (1991) found that the gas- 
phase CH 4 abundance is typically 10~ 3 times the CO abundance. If we assume the fractional 
abundance of CO to be 10~ 4 , this corresponds to the fractional abundance of 10 -7 . This 
value is lower than that observed in the present study by an order of magnitude. On the 
other hand, Boogert et al. (2004) reported the gas-phase CH4 abundance to be higher than 
3 x 10~ 6 toward the core of NGC7538 IRS9, and Knez et al. (2009) recently reported it to be 
4.8 x 10~ 6 toward NGC7538 IRS1. These abundances are comparable to our result, although 
their observation traces much hotter region (T rot =55-674 K) than our present observation. 

According to the chemical model calculation of dynamically evolving cores by Aikawa 
et al, (2008), CH 4 is mainly produced in grain mantles, and it is evaporated in a lukewarm 
region near the protostar. The gas-phase CH 4 abundance in L1527 comparable to or slightly 
lower than that expected from their chemical model calculation (10 -5 ). It should be noted 
that the fractional abundance of CH3D is also comparable to or slightly lower than the 
result of their new chemi cal model calculation involving the deuterated species (~ 3.6 x 



10~ 7 ) (lAikawa et al. 20121 ) . If the emitting region of the CH 3 D line is smaller than 20", the 



abundance of CH 3 D becomes higher, giving a better agreement with the model result. 

In this relation, it is interesting to compare our result with the solid CH4 abundance. The 
infrared absorption spectrum of solid CH4 is not available for L1527, because the protostar is 
heavily obscured even in the mid-infrared region. However, it i s detected in another WCCC 



source, IRAS15398-3359 flOberg et al. 20081 : ISakai et al. 2009af ). The abundance of the solid 



CH 4 relative to the H 2 ice is reported to be 0.06. When we assume the fractional abundance 
of the H2O ice relative to the gas-phase H2 to be 10 -4 , the fractional abundance of the solid 
CH 4 is estimated to be 6 x 10~ 6 . This is almost comparable to the gas-phase CH4 abundance 
found in LI 527. 
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Recent observational studies of low-mass star-forming regions reveal a significant chem- 
ical diversity among the low-mass Class protostars. Two distinct cases are hot corinos 
and the WCCC sources. The hot corino chemistry is characterized by rich saturated organic 
molecules like HCOOCH 3 and CH 3 CH 2 CN, which are evaporated from g rain mantles. The 



representative sources are IRAS 16293-2 422 and NGC1333 IRAS4A/B fe.g. lCazaux et al. 2003 



Bottinelli et al. 20041 ; ISakai et al. 20061 ). In these sources, carbon-chain molecules are gen- 



erally deficient. In contrast, the saturated organic molecules have not bee n detected so far i n 



the WCCC sources in spite of very sensitive single-dish observations (e.g. lSakai et al. 20071 ). 
Such a chemical variation originates from the difference in the evaporated species, and 
thereby in chemical composition of grain mantles. In the WCCC mechanism, the abun- 
dances of carbon-chain molecules critically depends on the amount of CH4 evaporated from 
grain mantles. If the gas-phase CH 4 abundance is higher than the abundance of OH (typ- 
ically 10~ 7 ; e.g. Herbst and Leung 1989), CH 4 can be a major destructor of C + instead of 
OH in a lukewarm region where H2O is still frozen out. Then, the WCCC occurs when the 
CH4 abundance exceeds this level. The CH4 abundance derived for L1527 certainly fulfills 
this requirement, which further strengthens the WCCC picture in L1527. 

In order to verify a role of CH 4 in the observed chemical diversity, it is important to 
measure the CH 4 abundances in various low-mass star-forming regions. As demonstrated in 
this study, observations of CH 3 D can be a useful method for this purpose. 
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Fig. 1. — Spectral line profile of CH 3 D ( = lo — Oo) observed toward L1527 with HHT. 
The total spectrum is the weighted sum of the V- and H-polarization spectra. The V- 
polarization spectrum shows a larger noise level than the H-polarization spectrum, because 
we have to exclude apparent bad-scan data caused by the receiver instability in integration of 
the V-polarization data. As a usual reduction process, each 5 minute integration scan were 
carefully checked to eliminate the bad scan data with a heavy baseline distortion and/or an 
apparent strong noise spikes. Note that no apparent spikes appear at the line channel. The 
on-source integration time for the H and V spectra are 25 and 13 hours, respectively. The 
vertical dotted line shows the averaged Vlsr value reported for deuterated species in this 
source (5.9 km s _1 ) (ISakai et al. 2009bl ). 
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Fig. 2. — Noise distributions of the spectra shown in Figure [U Panels (a), (b), and (c) 
represent the noise distributions for the total, H-polarization, and V- polarization spectra, 
respectively. This figure is prepared from the data for the central 93 channels excluding the 
line channels, where the baseline of the spectrum is subtracted. The distribution is close to 
the Gaussian form. The arrows indicate the peak intensities of the line. It is significantly 
detected for the total and H-polarization spectra. 
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Table 1: Parameters for the Observation of CH3D (Jk 



Parameter Value 

Telescope HHT 

7/MB 0.75 

HPBW 32 [arcsec] 

Rest Frequency" 232644.301 ± 0.075 [MHz] 

^mb 12 ± 2 [mK] 

/ T MB dv(3(r) 7.4 ± 2.8 [mK km s" 1 ] 

6.0 [km s- 1 ] 

rms 2.2 [mK] 



a The rest frequency is taken from CDMS (jMuller et al. 20051 ) . The uncertainty represents one standard 
deviation. 

b Obtained by a Gaussian fit. In this fit, dv is calculated to be (0.5 ± 0.1) km s" 1 , but it is affected by the 
backend resolution (0.32 km s _1 ). The uncertainty represents one standard deviation of the fit. 

c The rest frequency error (~ 0.1 km s _1 ) and the channel spacing of the backend (0.32 km s _1 ) are larger 
than the Gaussian fitting error. 
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Table 2: Column Densities and Fractional Abundances of CH 3 D and CH 4 . 



Temperature 


iV(CH 3 D) [cm- 2 ]" 


/(CH 3 D)*" 6 


iV(CH 4 ) [cm- 2 ] c 


/(CH 4 ) fe ' c 


15 K 


(6.3±2.4)xl0 15 


(2.1±0.8)xl0" 7 


(0.9-3.1) xlO 17 


(3.0-10.5) xlO" 6 


25 K 


(9.1±3.4)xl0 15 


(3.0±l.l)xl0~ 7 


(1.3-4.6) xlO 17 


(4.3-15.2) xl0~ 6 


35 K 


(12.9±4.9)xl0 15 


(4.3±1.6)xl0" 7 


(1.8-6.4) xlO 17 


(6.1-21.4) xHT 6 



a The error is evaluated from three times the standard deviation of the integrated intensities. 

b The values are calculated by assuming the H2 column density of 3 x 10 22 cm -2 . 

c The values for CH4 are calculated by assuming the CH3D/CH4 ratio of 2-7 % (See text). 



Note. - 



The source size is assumed to be 20" (jSakai et al. 20101 ) . 



